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Abstract 

The propagation of the scalar (a and 6) and vector {to and p) mesons in an iso-asymmetric 
nuclear matter is studied in detail, using the Walecka model. We calculate the invariant masses 
and spectral functions of the mesons, including the effect of meson mixing. At finite density, the 
mixing effect is quite important in the propagation of the scalar and (longitudinal) vector mesons. 
In the a channel, we find a three-peak structure in the spectral function, caused by the mixing 
effect. 
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The variation of hadron properties in hot and / or dense nuclear matter has recently drawn 
a great deal attention. In particular, the medium modification of the light vector (p, uj and 
0) meson masses, which may be a probe to study hadron properties in the extreme matter, 
has been investigated by many authors 

Experimentally, the photoproduction of uj mesons on a Nb target was recently measured 
at the ELSA tagged photon facility j^. The result shows that, for low momenta less than 
500 MeV, the in-medium uj meson mass at 0.6po (po is the saturation density of nuclear 
matter) is about 8% smaller than the free mass. Furthermore, the experiment to measure 
the invariant mass spectra of e"*"e~ pairs produced in 12 GeV proton-induced nuclear (C 
and Cu) reactions was also performed at the KEK proton synchrotron . On the low mass 
side of the uj meson peak, a significant enhancement over the known hadronic sources is 
observed. It is well understood by a model that takes account of the medium modification 
for the vector meson mass. These experimental results may suggest a downward shift of the 
vector meson mass in a nuclear matter (see also Ref. J]). 

On the other hand, the a meson has been treated as a correlated two-pion state in the 
scalar channel for a long time. However, some people have argued for the existence of the a 
as a genuine resonance 5| . Anyway, it is well recognized that the a meson is vital to produce 
the mid-range, attractive force between two nucleons in a matter. The 5 (or oq) meson also 
plays an important role in an iso-asymmetric nuclear matter. 

At finite nuclear density, the scalar meson can couple to the longitudinal (L) mode of 
the vector meson Gj. Furthermore, when charge symmetry is broken, the iso-scalar meson 
can couple to the iso-vector one. In an iso-asymmetric nuclear matter, where the number 
of protons, Z, is different from that of neutrons, N , the effect of charge symmetry breaking 
(CSB) is more enhanced than in vacuum, because the difference between the in- medium 
proton, Mp, and neutron, M*, masses is larger than that in vacuum [3, S]- 

This work is an extension and elaboration of studies in Ref. joj, where the iso-scalar a 
and UJ mesons moving in an iso-symmetric {Z = N) nuclear matter were considered. Here 
we generalize it in an asymmetric {Z ^ N) nuclear matter, including the iso-vector 6 and 
p mesons. We do not consider the propagation of pseudoscalar (tt and rj) mesons, because 
we use the relativistic Hartree approximation (RHA) in the calculation. Instead, we include 
the effect of the scalar or vector meson decay into these pseudoscalar mesons. 

We use a purely hadronic model, i.e., quantum hadrodynamics (QHD or the Walecka 
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model) [lo[. The Lagrangian density is given by 



(1) 



where Cqhd-i is the usual, QHD-1 Lagrangian density for the cr-cj-nucleon system [lOj and 
5CcT is a counter term for renormalizations. 

The 5 and p mesons are, respectively, described by Cs and Cp. The 5 meson couples 
to a nucleon through the isovector, scalar interaction: £™* = gsipTzSip, where ip (6) is the 
field of the nucleon (the neutral member of 6 meson) and gs is the 5-nucleon coupling 
constant. The nucleon mass in matter is thus given as MLn = Mfp\ — g^a =F gsS with a 
the a mean field and g^ the a-nucleon coupling constant. Here Mp{= 938.27 MeV) and 
M„(= 939.57 MeV) are, respectively, the proton and neutron masses in vacuum. The p 
meson has a tensor as well as a vector coupling to a nucleon, and the interaction is given 
by = —gpifj'^pP^Tz'ip + {fp/2Mi)%l)api,d'' p^r^tp , where i = proton (p) or neutron (n), pf^ 
is the neutral p meson field and gp (fp) is the vector (tensor) coupling constant. 

The pseudoscalar, tt and r], mesons are described by £7^7); in which their interactions to a 
nucleon are, respectively, written by the isovector and isoscalar, pseudoscalar couplings. 

Firstly, we have to solve the nuclear ground state within RHA. The total energy density 
is given by Stot = ^rha + '^^vf, where Srha is the usual one in RHA [lO|. The second term 
comes from the vacuum fluctuation correction (caused by the nucleon loop) to the a and 6 
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, ISj, and it is given by (also see around Eq.([7])) 
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where ma(s) is the a (5) meson mass, 
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1 — (4Mf/m^) with TTij the meson mass (j = a 



or (5), f{y) = 2v/=^tan-i ^-y-^ and M = (M^ + MO/2. 

The effective nucleon mass, M*, in matter is self-consistently calculated by solving the 
following, two equations: 
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where Ps{i) is the proton or neutron scalar density in matter and 
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A complete investigation of the propagation of a, cu, 5 and p mesons requires the inclusion 
of the meson mixing effect in a medium. To compute it, we sum over the ring diagrams ra, M, 

fin 

|8|, l9j, which consist of repeated insertions of the lowest order, one-loop proper polarization 
part. It is convenient to use the full meson propagator, Vab (a, & = 1 ~ 10), in the form 
of a 10 X 10 matrix. The lowest order meson propagator, is then given by a block- 
diagonal form: = block-diag (sq, Wq*^, c?o, Tq'^), where the free propagators for the a, u, S 
and p are, respectively, expressed by so{q) = Pail), Wo'^q) = ^^'^Puil), do{q) = ps{q) and 
''^o'^il) = i^^Vpil) with ^^^^ = —g'^'" + q^q^ lq\ and Pj{q) = {qf^ — m?j + ie)~^ (j = cr, c<j, 5 or 
p, and q^ is the four-momentum). 

Dyson's equation for V is also given in matrix form as 

v = v^ + p°np, (5) 

and the polarization insertion, 11, is written by a 10 x 10 matrix 

/ u^M) KM K&{q) KM) \ 
KAq) Kl{q) Ksi^i) KM) 
KM) KM K&{q) KM) 
\ KM) n^:^(g) ii%{q) KM) J 

where Ilafe (a, b = (J,uj, S, p) stands for the polarization for the process where the h meson is 
converted into the a meson through the nucleon loop. (Hereafter, Ilaa is simply denoted by 
H-a-) The polarization insertion, 11, can be separated into two subspaces, i.e., the iso-scalar 
and iso-vector block matrices. The coupling between the two blocks vanishes when charge 
symmetry is exact. However, even in vacuum, there exists the explicit CSB, i.e., the small 
difference between Mp and M„, and the CSB effect is magnified in an asymmetric matter j^. 



H 



(6) 



4 



Because the nucleon propagator, G(k), is divided into the Feynman (F) and density- 
dependent (D) pieces, the polarization due to the nucleon loop consists of the contribution 
depending on ps (the nuclear density) and that involving only the Feynman piece. The 
former (the D part) can be calculated analytically, while the latter (the F part) is divergent. 
We thus treat the F part using the method of dimensional regularization. For the detail of 

the D part, see Refs. Ban 

To regularize the nucleon loop contributions to the scalar meson propagators and the a-5 
mixing, we choose the counter term in Eq.([T]) as 

2 
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To fix the coefficients {Zj, Ai, Bi and Cim), we adopt the following renormalization condi- 
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where 11^ is the F part of the polarization insertion. This counter term yields the vacuum 
fiuctuation correction, Eq.([2]), to the energy density. 

In addition to the nucleon loop contribution, the lowest order polarization insertion for 
the a involves the pion loop contribution 9j 

n.,.(g) = -tgLml j ^A^(A;)A^(A; + g), (9) 

with g^T, the cr-vr coupling constant, the pion mass and ^-nik) the pion propagator. The 
imaginary part of Ila,-n describes the a meson decay into 2ti. In contrast, the 5 polarization 
insertion includes the r^-vr loop contribution [12|] 
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with qst^j^ the S-n-rj coupling constant and A^(A;) the r] meson propagator (m^ is its mass). 
To remove the divergences appearing in these polarizations, we use the renormalization 



conditions on the real parts of the polarizations: 
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For the vector meson, the F part of the polarization insertion arising from the nucleon 
loop, , again contains the divergent piece, and it may be removed by the condi- 
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The polarization insertion for the p meson also involves the pion loop 



[i2k + qY{2k + qYA^{k + q)A^{k) - 2^^'^A.(A;)] 



(13) 



where Qp^^ is the p-vr coupling constant. This again involves the divergence, but the same 
condition as in Eq. ffTT]) may remove it. 

For the p-uo mixing in vacuum, the polarization insertion due to the proton or neutron 
loop diverges. A finite mixing amplitude is then obtained by taking the difference between 
proton and neutron contributions p]. In contrast, the divergent piece arising from the 
tensor coupling is proportional to M* /Mi in matter, and hence it cannot be removed by the 
same method as in vacuum. However, we simply ignore this divergent piece in the present 
calculation.^ In the F part of the polarization insertion, the finite mixing amplitude from 
the tensor coupling is thus given by Ilpi*-'^'^ = {g^f p/'STi'^){Kp — K^), where 



M* 



X 



(14) 



The mixing between the scalar and vector mesons, e.g. Hau, etc., does not occur in 
vacuum and hence the regularization for such mixing is not necessary 8, 91. 

n Q 

The dielectric function, e, in a medium is now defined by [6|, 19| 



e = det(/ - P°n) = e| X esl, 



(15) 



^ The tensor interaction is not a renormalizable one in the conventional sense, and the strict renormaliz- 
ability is not necessary in effective field theories. The present prescription is equivalent to taking only 
the minimal, finite part of Ilpi*'''"'. This ambiguity may not affect much the final results of the invariant 
mass and the spectral function (see Figs. [U 2), because the ui branch is not close to the p branch in the 
dispersion relation and hence the effect of p-u mixing is small in matter. 
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where Etisl] is the dielectric function for the transverse (T) [scalar (S) and longitudinal (L)] 
mode. We find 



ST = (i- PpnJ)(i - p^ul) - Ppp.(njj2, (16) 
esL = [ii-pMi^-psT^5)-P.P5in^5f] [(i-ppnj)(i-p.n^;)-p,p^(n^j2] 

+ [(1 - pM{i - p^u^)pspMp)' + (1 - pMi'^ - Ppn^)P5P.(n°j' 

+ (1 - psUs)ii - p^n^Kpp(no^)2 + (1 - psUs)ii - ppn^^)p^pUKJ' 

+ 2p,p^nJ^|(i - pMpsIiLi^I + (1 - P5iis)p^ul^ul^} 

+ 2p^psU^s{ii - PpnJ')p^nLnL + (i - p.n^;)p,n°^n°^} 

+ 2p^psp,p^U^sT^liK.K + K^l) + ePMPpiK^K - K^l)'] ' (17) 
with ("^ = —q't/\q\'^, n"^(= = 11^^) the transverse polarization (we choose the direction 



of q as the 2;-axis) and 11^ (= 11'^'^ — n°°) the longitudinal one. The free propagator of the 
uj meson is now assumed to be = qfi ~ + ^^<^rJi, where we add the width of the 
uj in vacuum (Fj^ = 9.8 MeV) [9]. Using the dielectric function, one can calculate the full 
propagator, V (see Eq.([5])). 

In the present calculation, we use the following values of the free meson masses: = 983 
MeV, = 550 MeV, = 783 MeV, nip = 769 MeV, = 548 MeV and = 138 MeV. 
So as to fit the saturation condition of normal nuclear matter {£tot/ Pb — M = —16 MeV 
at po = 0.17 fm~'^), the symmetry energy [a^ = 30 MeV) and the observed p-u mixing 



amplitude in vacuum {{p\Hmix\^^)\qi=mi = —4520 MeV^) [l5|, we determine the coupling 
constants, g^, Qp and g^. The ratio of the tensor to vector coupling constants for the 
p meson is assumed to be fp/gp = 3.7, which is suggested by the vector dominance model. 
The coupling constants, gs^r] and gp^,, are, respectively, fixed so as to reproduce the decay 
widths of the 6 and p mesons in vacuum, = 59 MeV and 1° = 150 MeV. Furthermore, 
we assume the decay width of the a meson in vacuum, r° = 500 MeV, and determine the 
coupling constant, g^.,,, so as to yield it. We then obtain the following values of the coupling 
constants: g^^ = 8.119, g^ = 9.506, gp = 4.445, gs = 3.732, g^^,, = 23.65, gsnr^ = 0.4385 and 
Qp-K = 6.014. The present calculation gives the (proton-neutron averaged) effective nucleon 
mass, M* /M = 0.72, at po and the incompressibility, K = 465 MeV. 

Now we are in a position to show our results. The condition for determining the collective 
excitation spectrum is equivalent to searching for the zeros of the dielectric function. Because 
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FIG. 1: Invariant mass with or without the mixing effect (/p = 0.3 and q = 500 MeV). The left, 
two figures are for the meson masses without the mixing effect, while the right ones are for the 
masses including it. 

we are interested in the medium modification of the meson propagation, we restrict ourselves 
to the meson branch in the time-like region. 

In Fig. [H we present the "invariant mass" as a function oi ps- The mass is defined by 
m*„^ = \j — where go and q satisfy the dispersion relation for the meson branch. In 
the calculation, the proton fraction, fp = Z/{N + Z), is chosen to be 0.3, and the three- 
momentum transfer, q{= \q\), is 500 MeV. In the figure, we can clearly see the role of 
mixing. It is quite important in determining the meson mass at finite density. If the mixing 
is ignored, the masses of the L and a modes cross each other and the L mass of the p meson 
is below the other ones at high density. However, when the mixing is included, the L and 



8 



a masses never cross each other. At high density, the L mass is pushed upwards, while the 
a mass is pulled downwards. We can understand this phenomenon as a level-level repulsion 

n 

due to the mixing, which was first studied by Saito et al. [9]. It is familiar in conventional 
nuclear physics, e.g. in the Nilsson diagram. For the T mode, because the difference between 
the masses of the u and p mesons is large, the mixing effect is small even at high density. 
The similar result is obtained even in case of fp = 0.1. 

It is very interesting to calculate the spectral functions. Because, in a thermal model, the 
dilepton yields in heavy ion collisions are proportional to the spectral functions, the dilepton 
production rate per unit of four-momentum (d^q) is given by dNi+i- / d'^xd'^q = L^^H^'^, 
where d'^x is the space-time element, and L^y and H^'^ are respectively the lepton and 
the hadronic tensors in matter. The latter is then given in terms of the spectral function. 
The spectral function for the i-th mode {i = a, 6, uL, uT, pL, pT) is usually defined by 
Si{ml^v^ 9> Pb) = [Vi] /vr. 

In Fig. 2, we present the spectral function Si (at Pb/po = 2 and fp = 0.3) as a function 
of m*j^^ and q. The shape for the L or a mode is very complicated. In particular, the a 
mode is quite remarkable: at q=0 MeV there is only one peak, while at finite q three peaks 
appear in the spectral function. The two peaks, which have relatively large amplitudes, are 
seen at ~ m^j, and this two-peak structure is due to the a-u mixing [9j. In addition to 
it, a new, third peak appears around m*^^ ~ 0.4 GeV, which is caused by the mixing effect 
between the iso-scalar and iso-vector mesons. 

In the uL mode, the large peak is seen above m*^^ ~ 0.5 GeV, that approaches the usual 
w-meson branch in the limit ps 0. At finite density, the second, small peak appears 
below ^0.4 GeV. It is again generated by the mixing effect between the iso-scalar and 
iso-vector mesons, because such structure cannot be seen at fp = 0.5 For the pL mode, 
in addition to the main peak around m*^^ ~ 0.4 GeV, the second, small peak also appears 
above m*„^ ~ 0.5 GeV. In case of fp = 0.1, these second peaks become more clear. 

The spectral function is very simple for the T mode, because the mixing effect is small. 
However, even for the uT (pT) mode, the second, very small peak appears around m*^^ ~ 0.4 
GeV (above rn*^^ ~ 0.6 GeV), which cannot be seen in a symmetric matter j^. Thus, it 
is again due to the effect of iso-scalar and iso-vector mixing. In the 6 mode, the spectral 
function is very simple, because the mixing effect on it is quite small. 

The present result is not sensitive to the decay width of the a meson in vacuum. We 
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have checked it using r° = 600 MeV. In contrast, the spectral functions strongly depend on 
fp. In the (7 mode, one can clearly see that the third peak around m*j^ ~ 0.4 GeV grows 
very quickly as fp decreases. 

In summary, using QHD, we have studied the propagation of a, 6, uj and p mesons in an 
iso-asymmctric nuclear matter. We have illustrated that the effect of the meson mixing is 
quite important in the propagation of the a, longitudinal cu and p mesons at finite density. 
In particular, in the a channel, we have found a three-peak structure in the spectral function 
at finite three-momentum transfer. It is very interesting to measure such novel structure in 
future experiments. 
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